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1.  Introduction 


A  number  of  investigations  of  crack  problems  in  nonhomogeneous  media 
have  been  undertaken  in  which  the  elastic  moduli  vary  continuously  with 
spatial  coordinates-  In  all  of  these  studies  special  forms  of 
inhomogeneities  have  been  assumed  in  order  to  insure  a  tractable  problem 
for  which  the  asymptotic  form  of  the  stress  field  near  the  crack  tip  could 
be  calculated-  One  of  the  primary  objectives  of  these  works  was  to 
determine  the  effect  of  spatial  inhomogeneity  upon  the  known  singular  field 
quantities  associated  with  the  corresponding  homogeneous  problem. 

In  the  case  of  a  Mode  I  crack  embedded  in  a  nonhomogeneous  medium, 
symmetric  about  the  plane  of  the  crack,  specific  models  of  inhomogeneities 
examined  thus  far  show  the  usual  square  root  singular  crack  tip  stress 
field  associated  with  the  homogeneous  medium  occurs  for  the  nonhomogeneous 
problem  provided  the  shear  modulus  does  not  vanish  in  the  plane  of  the 
crack-  (see,  e.g.  Delale  and  Erdogan  (1983),  Gerasoulis  and  Srivastav 
(1980),  Schovanec  (1986),  and  Schovanec  and  Walton.) 

Studies  of  a  Mode  III  crack  located  on  the  interface  of  two  bonded 
materials  have  shown  that  the  assumption  of  spatial  inhomogeneity  in  either 
component  part  does  not  alter  the  known  square  root  singular  behavior 
present  in  the  piecewise  homogeneous  bimaterial  problem-  One  of  the  more 
general  models  utilized,  which  includes  the  homogeneous  and  bimaterial 

mediums  as  special  cases,  was  investigated  by  Delale  (1983)  who  took  the 

+ 

±  cx  x+g  y 

shear  modulus  in  the  form  p(x,y)  =  \x  e  with  the  crack  situated  in 

the  plane  y=0  and  U+,  3+(y  ,3  )  constants  specified  in  the  half  plane 
y>0  (y<0).  (See  also  Clements,  et  al-  (1978)  and  Dhaliwal  and  Singh 


(1978).) 
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In  the  situation  that  the  crack  tip  terminates  at  a  bimaterial 
interface  between  two  homogeneous  mediums  it  is  known  that  the  stress  field 


— Y 

singularity  is  of  the  form  r  ,  0  <  Y  <  1,  where  r  is  the  distance  from 
the  crack  tip*  (Erdogan  and  Cook  (1972))*  As  discussed  by  Atkinson  (1977) 

#  such  a  result  is  unsatisfactory  from  a  fracture  mechanics  point  of  view  in 
that  it  results  in  a  stress  intensity  factor  and  energy  release  rate  of 
zero  or  infinity.  To  alleviate  this  inadequacy  in  the  'ideal1  interface 

#  model  Atkinson  introduced  a  nonhomogeneous  layer  at  the  interface  in  which 
the  shear  modulus  varied  continuously,  matching  the  constant  values  at  the 
two  outer  mediums.  For  the  specific  forms  of  the  modulus  considered  it  was 


-1/2 

then  shown  that  the  stress  at  the  crack  tip  displayed  an  r 
singularity.  Assuming  that  the  elastic  module  could  be  expanded  in  a 
Taylor  series  about  r=0,  Atkinson  suggested  that  the  square  root  behavior 
would  always  result.  The  antiplane  shear  crack  perpendicular  to  an 
interface  was  also  considered  by  Erdogan  (1985).  To  investigate  the  effect 
of  a  continuous  but  nondif f erentiable  elastic  modulus  at  the  interface,  a 

3  x 

shear  modulus  of  the  form  y(x)  =  y^e  was  assumed  with  3  (3  ) 
corresponding  to  x>0  (x<0).  The  main  focus  of  that  study  was  to  show  that 
for  such  a  model  a  square  root  singularity  at  the  crack  tip  resulted. 
Erdogan  proposed  that  the  same  conclusion  should  be  valid  for  any 
continuous  but  nondif f erentiable  modulus  y(x). 

This  note  addresses  the  conjecture  raised  in  the  two  previous  studies. 


Without  assuming  any  particular  form  of  the  shear  modulus,  only  that  it  be 
continuous  throughout  the  medium  and  differentiable  everywhere  except  along 
a  curve  at  which  the  crack  tip  terminates,  it  is  shown  that  a  square  root 


3 


singular  stress  field  is  present  at  the  crack  tip.  As  will  be  clear  from 
the  subsequent  analysis,  the  situation  in  which  the  crack  lies  along  this 
interface  at  which  the  modulus  is  nondif f erentiable  is  handled  in  the  same 
manner,  so  that  in  such  cases  the  same  square  root  behavior  occurs.  This 
note  does  not  address  the  general  question  of  existence  of  solutions  for 
such  a  class  of  boundary  value  problems.  Rather,  the  approach  taken  here 
is  to  show  that  when  there  exists  a  physically  meaningful  solution,  a 
notion  to  be  made  precise  later,  then  that  solution  exhibits  a  square  root 
singular  stress  field. 

2.  The  Boundary  value  problem. 

The  specific  problem  considered  here  is  that  of  an  infinite 
nonhomogeneous  isotropic  elastic  solid  containing  a  semi-infinite  mode  III 
crack  terminating  along  a  curve  C  corresponding  to  the  interface  of  two 
bonded  materials.  Referred  to  cartesian  coordinates  x,y,z,  the  crack  is 
assumed  to  be  in  the  plane  y=0 ,  x<0  and  the  nonzero  component  of 
displacement  w(x,y)  is  related  to  the  nonzero  components  of  stress  by 
=  y3u)/3x,  0^^  =  p3u)/3y  where  the  shear  modulus  y(x,y)  is  assumed 
to  be  continuous  throughout  the  medium  and  differentiable  everywhere  but 
along  C.  The  primary  case  to  keep  in  mind  is  when  the  crack  is 
perpendicular  to  the  interface  in  which  case  C  corresponds  to  x=0.  If 
the  modulus  is  symmetric  about  the  plane  of  the  crack  then  the  particular 
boundary  value  problem  to  be  solved  is 


V 


div(y  Vo)  )  =  0 

✓  3(ji)  ,  _ x 

(x,0)  —  (x,0) 

3y 

(D(x,0)  =  0 


|x|<“,  y>0 
S(x)  x<0 

x>0 . 
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The  existence  of  a  solution  with  locally  finite  strain  energy,  i.e.  , 

J^(y/2)  |Vo)|^  dA  <  00  where  ft  is  any  bounded  measurable  subset  of  the 
upper  half-plane,  will  be  assumed.  Note  that  a  sufficient  condition  to 
insure  w  have  locally  finite  strain  energy  is  that  near  the  crack  tip, 

|Va)|  =  0(r  ^ ) ,  0<^Y  <  1.  The  point  here  is  to  show  Y  =  1/2.  A  growth 

condition  on  00  of  the  form  /  log 

|x  >R 


x  ( (Vp»Vo))/y  )(x)  dx  <  00  for  R>0 


is  also  assumed.  For  a  large  class  of  moduli,  (e.g.,  U(x)  =  (ctx+3)n  as 
considered  in  Atkinson  (1977))  this  condition  is  satisfied  if 

Vtu  =  0(r  ),  r*00,  which,  is  indeed  the  case  if  the  applied  tractions 

are  statically  self -equilibrating  (Cook  and  Erdogan  (1977)).  Again  it 
should  be  mentioned  that  the  question  of  determining  conditions  on  the 
applied  traction  S(x)  and  the  curve  C  which  assure  the  existence  of 
such  a  solution  u  is  not  taken  up  here. 

Define  the  Green1 s  function  by 


• 

r ( r , s ; x,y )  =  log  / [(r-x)Z  +  (s-y)Z]  [(r- 

-x)2  +  (s+y)  ] 

+  >Kr,s;x,y) 

where 

• 

— \  7  =0,  r  <°°,  s>0 

ar  as  1 

^(r,0)  =  log [ (r-x) 2  +  y2] 

r>0 

(1) 

• 

(r»°)  =  0 

r<0. 

(2) 

Applying  Green1 s  second  identity 
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9  v. 


v  V  2  v  V  j 

k  v1Av2  “  VVldA  =  ha  vl  9^  v2  TrT 


ds 


in  each  medium  determined  by  the  interface  C  with  v^  =  oj  ,  V2  =  F ,  and 
recalling  that  the  displacements  are  continuous  at  C  while  the  normal 
tractions  across  C  cancel,  there  results  the  representation  formula 


u>(x,y)  =  -/  /  r(s,r;x,y)  )  (r,s)  dsdv  +/  r(r,0;x,y)  dr* 


0 


s(r) 


V (r,0) 


With  the  above  assumptions  on  oa  the  first  integral  exists  and  the 
subsequent  expressions  to  be  derived  from  it  are  also  well  defined. 
Seek  the  harmonic  function  ^  in  the  form 


(3) 


’Kr,s;x,y)  =  ^  /  f(t;x,y) 


,  s2  ,  2 

(t-r)  +  s 


dt « 


Since 


||  (r ,0; s ,y)  =  I  /°°  ||-  (t;x,y)  ^  , 
—00 


inversion  of  this  relation  and  use  of  (2)  and  (1)  gives 


||  (r;x,y)  =  I  /“  ||  (t,0;x,y)  -°°<r<0 


t-r’ 


(4) 


and 


(r-x)  _  r“  9^  ,  n.„  ,  v  dt 

7 — 727T  Jn  (t’n’x,y)  t=F  * 

(r-x)  +y  0 


r>0. 


(5) 


Methods  for  solving  (5)  are  well  known  and  it  can  be  shown  that 


-1  f°°  ( u-x  )  / u  du 


o  w  —  \  r 

37  (t,0;x,y)  ‘  IZfl  ,  .2  2  u-t’ 

x  VT  0  (u-x)  +y 


t>0. 


(6) 
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Substitution  of  (6)  into  (4),  an  interchange  in  the  order  of  integration 
and  the  observation  that  for  u>0,  r<0 


leads  to 


r_L 


dt 


1 


✓  t  (u-t)<t-r)  /jTf11- 


|f  (r;x,y)  -  f  ^  <U'X)  d“ 


2/- — -  i.  ,  v2,  2  u-r  * 

7T  /  |  r  |  0  (u-x)  +y 


r<0 . 


An  integration  by  parts  gives 


14  =  1 


r 


4?r  /  |r  |  0  / u  (u-r)' 


fr-^  log[(u-x)2+y2]  du, 


r<0.  (7) 


At  this  point  it  is  convenient  to  observe  that  if  (7)  is  differentiated 
with  respect  to  y  and  then  integrated  with  respect  to  r,  keeping  in 
mind  that  f(* ,x,y)  must  vanish  at  infinity,  then  one  has  for  r<0. 

(  l  l  /IFF 


9f 

97  = 


IT  / 


x  x-r 


0 


x>0 


x<0 


(B) 


In  order  to  illustrate  the  singular  nature  of  the  stress  field  and  to 
determine  the  stress  intensity  factor,  only  3o)/8y  (x,0)  is  required.  The 
method  of  computing  is  carried  out  in  the  same  manner. 

3.  The  stress  singularity  and  the  SIF. 

From  the  definition  of  T,  (1),  and  the  representation  of  o>  given 
by  (3),  one  has  that 


(x,0)  =  lira  {J  J  {J  1L  (t;x,y) 


3oo 

77  ^x’u;  =  11™+  u..  j  u  77 


y+0  -°°  0 


(t-r)  +s 


2 — 2  dt^  (7y^7ai)(r,s)  dsdr 


(9) 


7 


+ 


1/ 

IT  J 


y 

2  2 

(r-x)  +y 


s(r) 

U(r,0) 


0 

dr  +  / 


(r;x,y) 


s(r) 

M(r,0) 


dr}  . 


Making  use  of  (8)  it  easily  follows  that  the  limit  of  the  single  integrals 
in  (9)  is  given  by 


'  s(x) 
P(x,0) 


i 


l71 


i  0 

Yx  “°° 


s(r) 

u(r) 


dr 

r-x 


x<0 


x>0. 


(10) 


By  utilizing  (8)  and  (1)  it  can  be  shown  in  a  similar  manner  that  the  limit 
of  the  triple  integral  in  (9)  is 


j-  r  n 


7T  '  '  1  ,  .2  2  7T  /  — 

-°°  0  (x-r)  +s  /x 


+  I  I  rO  /jtl 


x-t  ,  v  2  2 

(t-r)  +s 


dt} 


Vp*  Vco 


dsdr,  x>0 


(11) 

x<0 . 


Combining  (10)  and  (11)  it  follows  that  for  x>0, 


Vjj  *  Vo) 


o  (x,o) .  r  r _ v-t— 

’  -0  <*-r)2+s2  " 


dsdr  + 


aiziSl  r  1“  1°  ^  — £_  dtdsdr  +  HOiiill  / o  /T7  s<r2 

Tr/x  0  X_t  (t-r)2+s2  y  W-  _co  1  V(v,0 


X 


(12) 

dr 
0)  r-x 


•  From  (12)  the  square  root  singular  behavior  is  apparent  and  in  the  case 

that  p  is  constant,  the  known  result  for  the  homogeneous  medium  (Willis 
(1967)), 

ayz(*>°)  =  /°  /TH"  s(r) 

yz  t i/x  -°°  r  x 

is  easily  recovered.  Since  a  square  root  singularity  is  always  present  at 


the  crack  tip,  the  stress  intensity  factor  may  be  defined  as 
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K  =  lira  y(x,0) 
x+0 


✓5 


3o) 

■37 


(x,0) 


and  from  (12)  one  has 


k  =  "Jp  (  1°  —  uTToT  dr  +  r  /“/”  — - Va  U 

,0>  -*  0  /“[Tf  (t-r)2+s2  “ 


Vp*Vcj  1 

— - -  dtdsdr  j  • 


As  pointed  out  above,  the  calculation  of  O ^  is  carried  out  in  a 


similar  manner  and  by  changing  to  polar  coordinates  it  is  not  difficult  to 


show  that  in  general  o..(r,0)  =  o(  — ) ,  r+0.  Finally  it  should  be  noted 

ij  /r 

that  a  similar  derivation  can  be  carried  out  for  finite  domains.  In  this 
case  the  appropriate  boundary  conditions  must  be  prescribed  and  the 
analogous  Green’s  function  computed.  The  growth  condition  imposed  on  oa 
however  would  no  longer  be  necessary. 
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